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t****" ' The Bogoliubov equations of the quasi-particle excitations in a weakly interacting trapped Bose- 

OS . condensate are solved in the WKB approximation in an isotropic harmonic trap, determining the 

Os ■ discrete quasi-particle energies and wave functions by torus (Bohr-Sommerfeld) quantization of the 

integrable classical quasi-particle dynamics. The results are used to calculate the position and 
strengths of the peaks in the dynamic structure function which can be observed by off-resonance 
O ' inelastic light-scattering. 



. I. INTRODUCTION 

The experimental realization of Bose-Einstein condensates of weakly interactingalkali atoms confined by magnetic 
traps has revived the interest in the properties of Bose-condensed systems || , in particular in solutions of the 
Gross-Pitaevski ^ equation in an external potential, describing the condensate, and of the Bogoliubov equations 
governing the quasi-particle excitations out of the spatially inhomogeneous condensate . Approximate spatially 
inhomogeneous solutions of the Gross-Pitaevski equation have been obtained by a Gaussian variational ansatz, which is 
appropriate for moderate size condensates, and by the Thomas- Fermi approximation G3|, which is very good for large 
' condensates, except in a narrow region around the surface of the condensate. Similarly, solutions of the Bogoliubov 
equations for the low- lying quasi-particle excitations have been constructed numerically ffi], by a variational method 
Q , and in the hydrodynamic limit || . Such low- lying states have recently been excited by resonant periodic variations 
in the trap parameters jll] --|l3), and good agreement with the theoretical predictions has been obtained. 

Much less work has been done so far on the high lying excitations, where all the methods which have so far been 
used for the low lying ones are either impracticable or don't apply. Such high lying states are of experimental interest 
because they are excited by off-resonance inelastic light scattering. In a previous paper |fl4|| we have computed the 
high- lying quasi-particle energies and the amplitude of the inelastic peaks in the scattering function using the Thomas- 
Fermi approximation and solving the Hartree-Fock-Bogoliubov equations in perturbation theory of first order for the 
energies (with an estimate of the second order ones) but only to zero order in the wave functions. It is desirable to 
improve on this highly simplified approach by developing a self-consistent non-perturbative approach and extend the 
calculation to the whole energy region. It is our purpose, in the present paper, to present a non-perturbative method 
— the WKB approximation — for the solution of the Bogoliubov equations in a trap which allows us to interpolate 
between the hydrodynamical and the high-energy regime, and to use it to calculate the energies and wave-functions 
of the quasi-particle states. In a previous brief report we have already presented some results of our calculation of 
the high lying levels of the quasi-particles fLq] . Here we shall apply our results to calculate the dynamical structure 
function, which is of high interest, because it could be measured in off-resonance light scattering p6|-p8|. Calculations 
of the dynamical structure function have previously been presented for the homogeneous ideal Bose gas Jl8|,|l9| and 
for the homogeneous interacting Bose gas ^0|, which are however not directly applicable to the trapped condensates. 
Results have also been presented for the hydrodynamic limit of the structure function [ plf . Very recently results 
for the dynamical structure function in the local density approximation have been presented ^2j, However, due to 
the limitations of that approximation, the resulting structure function cannot show structure beyond the constraints 
imposed by the locally averaged energy and momentum conservation. 

The results we shall present here are applicable in a large enery domain ranging from low lying to high lying levels 
of the trapped condensates and retain the full information on the discrete nature of the quasi-particle levels. They 
are, however, restricted to the spatially isotropic case. In previous work |l4| we have calculated the dynamical structure 
function for the very high lying quasiparticle levels in isotropic trapped condensates in a simpler approximation, 
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replacing the high lying states by free trap states. The more sophisticated non-perturbative calculation presented in 
the present paper shows, however, that the interaction of the high lying states with the condensate has a non-negligible 
influence on the selection of the angular momentum of the states excited by light-scattering. 

An essential and independent input for our WKB-calculations is the condensate wave-function, which must be 
obtained by solving the Gross-Pitaevski equation. As we intend to carry out the WKB calculations analytically, as far 
as possible, we shall use here the main analytical approximation for large Bose-condensates with repulsive interaction, 
the Thomas-Fermi approximation. The necessary results are collected in section |J. 

In section III. the WKB approximation for the Boguliubov equations in an isotropic harmonic trap is developed. It 
is implemented in section |Tvj . by studying the classical quasi-particle dynamics, which is integrable in isotropic traps. 
For a discussion of the non- integrable quasi-particle dynamics in anisotropic traps see p3| . 

The reduced action is calculated as a function of energy and angular momentum. In section |Tv| . the Bohr-Sommerfeld 
quantization of the energy levels of the classical quasi-particle dynamics is performed and the semiclassical wave 
functions, properly normalized, are obtained in section |v|. The necessary integrals, though tedious, can all be done 
analytically. 

In section |vi| . we present our results for the scattering functi on and in section VI I some numerical examples are 
given. Our conclusions are summarized in the final section VIII . 



II. BASIC EQUATIONS FOR CONDENSATE AND QUASI-PARTICLE EXCITATIONS 

A necessary prerequisite of any calculation of the quasi-particle excitations is an appropriate description of the 
condensate. The relevant results are summarized in the present section. For weakly interacting atoms at zero 
temperature practically all particles are in the condensate and the energy £ is a functional of the macroscopic 
condensate wave function ipo( r ) 



— \VMr)\ 2 + U(r)\Mr)? 



2nh 2 c 



in 



\Mr)\ 4 



(2.1) 



where the constraint 



d 3 r\Mr)\ 2 = N 



(2.2) 



must be imposed. U (r) is the trap potential. In our calculations we shall assume an isotropic harmonic trap 



U(r) = ^mu^r 2 



(2.3) 



for simplicity. The scattering length is assumed to be positive and is denoted by a. The number of atoms in the 
condensate is denoted by No. We assume that the interaction is very weak and the temperature is close to zero, so 
that Nq ~ N. The mi nim ization condition for E, introducing the chemical potential \x as a Lagrange multiplier to 
enforce the constraint (2.2), yields the Gross-Pitaevski equation || 



— V 2 Vo + U{r)ip a + 



(2.4) 



A. Thomas-Fermi approximation 



In the present paper we shall use as an analytical approximation for the minimum of (2.1) the Thomas Fermi 
approximation, which holds in the limit of large condensates Nqcl » do = 



good approximation to solve (2.4) by neglecting the kinetic energy term |K 



Mr) 



1/2 

Qi-U(r))) 9(/i-C/(r)) 

47m a 



In that limit it becomes a very 



(2.5) 



where \i must be determined as a function of No by imposing the normalization condition (2.2), which yields 
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M=^(15iY aM)) 2/5 (2.6) 



The conditions [i » Tiujq and Nqcl >> da are therefore equivalent. In eq. (2.5) 0(x) is the step function. The 



approximation eq. ( |2.5| ) with discontinuous derivative at r = r = y^/i/rowf is sufficient, except in those cases where 
derivatives of the wave- function ipo( r ) for r close to ro are important. In the latter case a boundary layer theory can 
be used to compute the condensate wave function close to the Thomas Fermi radius , but we shall not need this 
improvement in the present paper. 

B. Quasi-particle excitations 

Let us now turn to a Bose condensate with excitations and adopt a description in the Heisenberg picture, where 
the field operator satisfies the equation of motion 

^d'tb Ti 2 7 ,_. . T 4:TTh 2 a 7,77 .„ 
at 2m m 

For a weakly interacting Bose condensed gas the Bogoliubov approximation consists in splitting the total Bose field 
operator ip(t) into 

m = fofo + ^(t))e-^l h , (2.8) 

where ip is the time-independent condensate wave function and <p(t) describes the excitations from the condensate, 
and to linearize the field equation for 6{t) 

lh ^T = " ^ V + U r - M + hAo r 2 + ifoto V . 2.9 

at \ 2m m I m 

The Bogoliubov transformation 

= E' (W^e"^* - v*(r)ap^A (2.10) 
j 

and its adjoint, with Bose operators a>j, for the quasi-particle excitations and with the sum over j restricted to 
states with ojj 7^ 0, yields the Bogoliubov equations 



-|^V 2 + U cS (r) - E j -K(r) 

-K*(r) -^ 2 + U cS (r)+E J 



2m 

where we used the abbreviations 



(2.11) 



8wh 2 a, 



Ej = huj , U eS (r) = U{r) + — — (r) \ A - \i 
K(r) = ^Mr) 2 - (2.12) 



Eq. ( 2.11 ) is consistent with and its solutions must be chosen to satisfy the orthonormality conditions 

d 3 r(u jU * k - VjV* k ) = S jk 

d 3 r(u*v k - u k v 3 ) = (2.13) 



in order to guarantee the Bose commutation relations of the ay, aQ~. 

A formal solution of eq. ( 2.11 ) at zero energy Ej = is given by the condensate 



i,=v*=Mr) , Ej=0, (2.14) 
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but this solution is not normalizable in the required sense. For a harmonic trap U (r) = Y^=i ^o^l three linearly 
independent and correctly normalized exact solutions 

Ui = bfip (r) , Vi = biii* {r) , E { = hu> ai 

always exist pq| , where 

h d 



(2.15) 



h = 



mujQi 



2mu)oi dxi V 2Ti 
are the Bose operators for the harmonical oscillations in the trap 



Xi , bj = (bi)~ 



(2.16) 



III. WKB-APPROXIMATION FOR QUASI-PARTICLE STATES IN THE CLASSICALLY ALLOWED 

REGION 



We shall here treat the WKB approximation for the case of an isotropic harmonic trap, where U e g(r) and K(r) 
depend only on the radial coordinate r. 

Let us drop the index j for the time being and expand formally in h 



u = (a„(r) + h ai (r) + h 2 a 2 (r) + . . .)e iS ^' h Y lrn {e, 4>) 
v = (6o(r) + hh(r) + h 2 b 2 (r) + . . .)e iS(r)/fi ^ m (0, 0) 



(3.1) 



In Eq. (3.1) we shall treat U e g(r) and K(r) as being independent of fi. The angular dependence is separated out by 
the spherical harmonics Yi m . The angular momentum, denoted by J in the following, is then quantized in the ususal 
way by J 2 = Ti 2 l(l + 1), w hich will b e replaced in our WKB-treatment according to Langer's rule J 2 = fi 2 (l + 1/2) 2 . 
Inserting the ansatz (3.1) into (2.11) we obtain up to order Ti 2 the equations 



with 
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-K(r) 
$- + tAt + U eS (r) + E 
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where p r = dS/dr. 



(3.2) 
(3.3) 
(3.4) 

(3.5) 
(3.6) 
(3.7) 



The solvability condition of (3.2), the vanishing of the determinant of L , gives the condition 



det(i ) 



2m \ dr 



J 



2mr' 



\K(r)f 



0. 



(3.8) 



We can solve this equation for E, restricting ourselves to the branch E > 0, because the replacement E — > — E merely 
amounts to the renaming u — > v* , v — > u* . Then we obtain the time-independent Hamilton- Jacobi equation 



H(dS/dr, r)=E 



(3.9) 



with 
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H(Pr,r) = \ £- 



I'T 



J 2 



+ U eS (r) -\K(r)\ s 



(3.10) 



2m 2mr 2 

This defines the Hamiltonian for the classical quasi-particle dynamics. Introducing the time-dependent action 

W(r, t) — S(r) - Et (3.11) 
we may replace eq. ([T|) by the time-dependent Hamilton Jacobi equation 



dW TT ,dW . 



o. 



(3.12) 



It follows from eq. (3.8) that with S also —S is a solution of the Hamilton Jacobi equation ( p.9[ ). It is interesting to 
note that the classical quasi-particle dynamics is integrable only in the special case of an isotropic harmonic trap. Even 
in the experimentally realized case of an anisotropic axially symmetric trap the classical quasi-particle dynamics turns 
out to be nonintegrable, in general. A detailed investigation is given in P3fl . In this case a semi-classical quantization 
of the quasi-particle dynamics requires the methods introduced in the field of quantum chaos, such as semiclassical 
quantization based on periodic orbits, or even replacing the quasi-particle spectrum and wave-functions by results 
based on random matr ix t heory. 



After satisfying eq. (3.8) the general solution of eq. (3^2) can be written in the form 

(J) = ao(r) ( 1 z E±V^EMl) 



(3.13) 



where ao(r) is still arbitrary 

Lq is not invertible. The solvability condition of eq. ( |3.3|) is that its inhomogeneity is orthogonal to the kernel of 
L 



{a* ,b* )Li 



0. 



(3.14) 



Evaluating this condition for 00(7") ^ we obtain after some rearrangements a classical transport equation for ao, 



dS_d_ 

dr dr 



In 



(|ao(r)| 2 (l + Wl + E 2 /\K{r)\ 2 -E/\K{r)\f 



1 d 

r 2 dr 



which can be rewritten as a continuity equation 

r 2 \a (r)\ 2 (l + (^1 + E 2 /\K(r)\ 2 ~ E/\K(r)\) 



_L_ d_ 

r 2 dr 



2\ dS_ 
dr 



,dS 
dr 



0. 



= 0, 



(3.15) 



(3.16) 



The physical meaning of this result as classical conservation law for the quasi-particle current becomes clear if we 
rewrite it once more as 



l__d 

r 2 dr 



[r 2 (|ao(r)| 2 - |& (r)| 2 H(r-)] = 0, 



where we used the identities 
d 



v r {r) 



dp r 



H(p r ,r] 



Pr =dS/dr 



^E 2 + \K(r)\ 2 1 dS _ 1 \K\ 2 + {^E 2 + \K\ 2 ~E) 2 dS 
E mdr ~ m \K\ 2 - (^E 2 + \K\ 2 - E) 2 dr 



(3.17) 



(3.18) 



to e xpres s dS/dr in terms of the radial component v r (r) of the quasi-particle velocity field. The general solution of 
eq. ( 3.17| ) can therefore be written as 



|a | 2 H&o| = 



const 
r 2 v r (r) 



(3.19) 



In principle, this procedure can be continued. Having satisfied the solvability condition for eq. (3.3) we can determine 
its general solution in the form 
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= oi(r) 



-fi/l^l + vT+W 



( & f rt (r; 



(3.20) 



Here the first term is the general solution of the hom ogeneous equation whose still arbitrary amplitude a\ (r) must be 
determined from the solvability condition of eq. (3.4) 



(a* ,b*){L 



«i 



= 0. 



(3.21) 



The second term in eq. (3.2C) is a particular solution of eq. (3.3), which we choose with vanishing amplitude a^ art (r) = 
in order to define the first term in eq. ( 3. 20]) in an unambiguous way. As we shall not need the explicit form of oi, b\ 
we shall here not work it out from eqs. (3.2C), (3.21). 



IV. SEMICLASSICAL QUASI-PARTICLE ENERGIES IN ISOTROPIC HARMONIC TRAPS 

For the radial canonical momentum as a function of r we obtain from ( 3.1 0| ) 



p r = V2m\ y / E 2 +K 2 (r) - 



J 2 
2mr* 



- U cS {r) 



(4.1) 



The semi-classical quantization can now be performed in a straightforward manner by putting for the angular mo- 
mentum 



and for the radial action variable 



I r (E, J) 



J = Je = h(£+l/2) 



p r (E, J,r)dr = h(n r + 1/2) 



(4.2) 



(4.3) 



Here £, n r are the integer angular momentum and radial quantum number, respectively. W e us ed 'Langer's rule' p6| of 
semi-classical angular momentum quantization by including 1/2 on the right hand side of (4.2), and took into account 
a phase-shift of 7r/4 in the wa ve function at the classical turning points r<, r > of the radial motion by including 1/2 
on the right hand side of (4.3). 



For the Thomas Fermi approximation of the condensate the radial integral in eq. (4.3) can be performed analytically, 
even though the resulting expressions become rather lengthy. In this manner an implicit analytical formula for the 
energy levels in WKB approximation is obtained. Numerically solving the implicit expression for E for different values 
of n r , £ we obtain the energy levels E rir i in WKB approximation. 



A. Types of motion 

In Thomas-Fermi approximation (|2.5|) we have 



1 2 2 

-uj n r 



U eS (r) 

X 2 (r) = (M-f^)^e( M -^ V) 



(4.4) 



where 



hw 2 
— ^-Kn 



Rn — 



^ 15jV oa y /5 r 



do 



mujQ 



(4.5) 



The Thomas- Fermi approximation requires Rq >> 1 for consistency. Let us introduce scaled variables and parameters 
by 



G 



p=(r/d R ) , e = E/p , j = (J/fiflg) , 7T r = Pr (d /hR ) . (4.6) 



Then eq. ( fl.l| ) becomes 



^ e 2 + (1 _ /9 2 ) 2 8(1 _ p) _ pi p 2 _ |! _ p 2| , (4 . 7) 



In order to find the turning points p<, p> we have to look for the solutions of 7r r = 0. 
i) Domain p < 1: Turning points in this region have to satisfy 



7 2 

2 J 
P< = ^ 



2.r 



(l + Vl + e' + Sj 2 ) . (4.8) 



The negative root is dismissed because only real and positive solutions are physical. The remaining solution (4.S) 
turns out to be a lower boundary of the p- values for which 7r r is real. The condition 

p\ < 1 (type B motion) (4.9) 

is only fulfilled if 

e > j 2 (tyP e B motion) . (4-10) 

We shall denote this motion as being of type B. As there is only a lower turning point p < within the region p < 1, 
the upper turning point must lie in the region p > 1, i.e. for type B motion the quasi-particle leaves and reenters the 
condensate during each radial period. 

ii) Domain p > 1: Turning points in this region must satisfy 



P> 



1, _ 1 



j ±(e + l)±-V(e+l) 2 -4f . (4.11) 
They exist if and only if 

e > 2j - 1 (4.12) 



where we recall that e is restricted to e > by definition. The lower turning point given by eq. (4.11) exists only if it 
satisfies 

p\ > 1 (type A motion) (4-13) 

in which case the classical quasi-particle dynamics is entirely confined to the region outside the condensate. It is 
then simply the motion in the harmonic potential of the trap. We call this motion being of type A. It requires the 
conditions 

e > 1 , j 2 > e > 2 j - 1 (type A motion) (4-14) 

all to be satisfied, see fig. [l]. 

If the turning point p< does not exist in the region p > 1, then the turning point p < must occur for p < 1, i.e. the 
only alternative of type A motion is motion of type B satisfying eqs. ([1.9]), (4.1C). The upper turning point of both 
types of motion lies in the region p > 1 and is given by 



,2 



-(e+1) + -\/(e+ l) 2 - 4j 2 (type A, B motion) . (4.15) 



It is easily checked that p^ > 1 for both types of motion. In fig. [j] the parts of the angular momentum — energy 
plane accessible to the two types of motion are shown. Let us now proceed to evaluate the radial action of both types 
of motion and quantize. 
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B. Asymptotic limit of a large condensate Rq — > oo 



The integral to be evaluated is 

fiR 2 r p> 

' n r dp. (4.16) 



In order to obtain simple closed expressions we first study the case of a large condensate for which R — > oo, while 
we keep the energy E and angular momentum J of the quasi-particle fixed, i.e. the scaled quantities 



^=—152 - 3 = —of- ( 4 - 17 ) 



2E 1 . _^ + l/2 

are both 0(i?Q~ 2 ) in that limit. It follows that the necessary condition for type A motion e < j 2 (i.e. motion which 
remains outside the condensate) cannot be fulfilled in this limit, and therefore only type B motion can occur. In the 
present limit we have asymptotically 

P< ^ , P> ^ Vl + e (4.18) 

^/2E 2 + (hu; ) 2 (£ + 1/2) 2 



+e(p- i)^ + i-p 2 

The integral 



hR_f 

7T 

can be performed asymptotically for e — > with the result 



I > = —^ y/l + e- p 2 dp (4.20) 

7T ./i 



tl i f2E^ 3/2 



The contribution of the part of the type B motion outside the condensate to the total radial action therefore becomes 
negligible 0(Rq 1 ) as Rq — ► oo. The 0(l)-part of the action is therefore given by 



3* i .£±11!,,. (4 . 22) 

7T 7 P< \ hwl (i - p 2 ) p 2 



The substitution 



p2 = l - ww^IW^Wn sin2 v (4 ' 23) 



reduces the integral to a rational function of sin <p and we obtain after integration 

h 



Ir 

2 



^(E/hw ) 2 + (i + l/2) 2 -(£+ 1/2)) + (J^j . (4.24) 

From the semi-classical quantization rule I r — h(n r + 1/2) we obtain the energy levels 

E nrt = huj (2n 2 r + 2n r l + 3n r + i + l) 1 ' 2 + ^ . (4.25) 

This result is remarkably close to the hydrodynamic result 

£;hyd _ ^ ( 2n 2 + 2n r £ + 3n r + £) 1/2 (4.26) 

derived by Stringari || for the low lying levels in the hydrodynamic regime. In that regime the WKB-method is 
not applicable, but for energies E Hr £ >> Tiloq the agreement between (4.25) and (4.26) becomes very good. Thus 
we find that there is a common regime of applicability of both the hydrodynamic approximation and the WKB 
approximation, where the energy is sufficiently large to apply the WKB approximation but not too large to invalidate 
the hydrodynamic approach. 
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C. Large-energy limit 



For type A motion the classical particle never enters the condensate and the Bohr-quantization (4.3) leads to 

E rirl = hiu (2n r +£ + 3/2) - fi = E° sc . (4.27) 



which depends only on the pricipal quantum number n — 2n r + 1 (see equation ( 4.35| ) below). However, for type B 
motion, i.e., e > j 2 , E nr £ differs from E° sc e 



E n 



El 



(4.28) 



The shift 8 n i lifting the degenerecy of the free oscillator levels is expected to be small in the large-energy limit 
En r e ^ p and can be calculated in the following way: 
For type B motion (|4.3|) can be written as 



1/2 = ^ f P> n r (p) dp=^ [ dpyj V^ 2 + (1 - - f/p 2 ~ (1 - V) 



ri r 



^ /'' dp^ ( ^3 2 /P 2 + ^-P 2 )-— I We 2 -J 2 /P 2 + (1~P 2 ), 



2 r 1 



7T J p M 



IT J p (°) 



(4.29) 



where /?< and p > are the turning points for 7r(p) = \/e 2 — j 2 /p 2 + (1 — p 2 ), i.e., for the free oscillator: 



(o) 

P> ,P< 



\ 



e + 1 



± 



e+1 



(4.30) 



The second integral in ( 4.29| ) is trivial and using (4.28) eq. (4.29) can be written as 

2R 2 n 



2R 2 



2 r 1 



d P ^ V* 2 + (i - f/p 2 - (i - P 2 ) 



dpy/e 2 -j 2 /p 2 + (l-p 2 ). 



(4.31) 



For type B motion j 2 < e, thus in the present limit e — > oo and a = j 2 /e is kept fixed (0 < a < 1). Our small 
parameter is 1/e. In this limit the two integrals in leading order cancel each other. With the replacement p 2 = 1 — z 
the next to leading approximation gives 



R 



l-o 



dz 



VI — z — a iiTyfl 



4jR l(l_ a) 3/2. 



(4.32) 



The shift is clearly tending to zero in the large energy limit, thus one can replace the unsealed energy E nr £ with E° s 
on the right hand side of eq. (4.32): 



5 n (£) 



1 [(4 M /^ )(n + 3/2 - iJ,/huj ) - (£ + 1/2) 2 )] 
3tt [ n + 3 /2_ n/hujo] 2 



3/2 



(4.33) 



which is the asymptotics of the shift in the large energylimit. This is the result presented in [|15|. For large values of 
n = 2n r + 21 it agrees with the perturbative result of |L4[] . 



D. Radial action in the case of a finite condensate 



We now evaluate the radial action for finite Rq- The two types of motion must be considered separately. 
Type A motion: 

This type of motion requires 1 < 2j — 1 < e < j 2 . Its radial action is that of the free harmonic trap and easily 
obtained as 
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Ir(e,j) 



fiRl f p > dp 



7T J p< P 



(4.34) 



where p > , p < are given by the two roots (4.11). Quantizing according to eqs. (4.2), ( |4.3| ) and solving for the energy 
levels we obtain 



E„ rt = hu (2n r +£ + 3/2) - p 



(4.35) 



i.e. the result which is obtained for the free trap, shifted by the change of the ground state energy due to the 

condensate. 

Type B motion: 

For this type of motion we have j 2 < e. The radial action is given by the integral 



Ir(e,j) 



Mil [P> 



dp ( e 2 + (i- P 2 ) 2 e(i-p)) 



,1/2 _ 3_ 

' P 2 



I1-P 2 I 



1/2 



(4.36) 



where p <: p > are given by eqs. (4.8) and (4.15), respectively. 

The integral can be performed in terms of elementary functions. We give a summary of the necessary steps in the 
appendix. The result is most succinctly expressed in terms of the functions 



^ 2 COS 2 if 



a + b sin tp 



dip 



J> (A, a, b) a 2 > b 2 
J< (A, a, b) a 2 < b 2 



(4.37) 



^ 2 co S 2 p f K> (A, a, b) a 2 > b 2 

A (a + 6sin^^-\^<(A,a,6) a 2 < b 2 ^ S8 > 

which can be written and programmed explicitely in terms of elementary functions. The radial action is then obtained 
in the form 

HR 2 f 3 

I r (e,j) = ^2c i J > {A 1 ,a i ,b 1 ) + C 4 J<(A X! a 4 , &i) 



with the parameters 



. i=l 



+C 5 K < (A 1 ,a i ,b 1 ) + C 6 J>(A2,a 5 ,b 2 ) > (4.39) 



Ai = arcsin [ — , 6 + J ' | (4.40) 
A 2 = arcsin f 1 ~ C = I (4.41) 



v /( e + l)2_ 4j 2 



61 = ev 7 ! + e 2 + 2j 2 (4.42) 

6 2 = Jf^V-i 2 (4-43) 



ai, 2 = T(e + J 2 )v / e 2 + l+/(e-l) (4.44) 

a 3 = -e - e 2 - 2j 2 (4.45) 

a 4 = e - e 2 (4.46) 

as = ^ (4.47) 



and the coefficients 
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C\.2 

c 3 



I(l-eTV^TT) 

le 2 (l + e 2 + 2j 2 ) 
4 



e(l + e 2 + 2j 2 ) 



1 



e(l + e 2 + 2j 2 ) 



^4=^(6- 2) „ 

1 



C 5 = -^TjV(l + e 2 + 2j 2 ) 



(4.48) 
(4.49) 
(4.50) 
(4.51) 



(4.52) 

The semiclassical approximations for the quasi-particle excitation energies are now obtained by solving the equation 

/ 2E l + l/2\ , . . 

n^'^r) =7lK+1/2) - ( 4 - 53 ) 

for E. This last step cannot be done analytically, but it is easily accomplished numerically. 



V. SEMICLASSICAL WAVE FUNCTIONS 



For purely radial dynamics eq. (3.1£) reduces to 



1 d r 

r 2 dr 



2 2 (r)r 2 (l + {^l + E 2 /K 2 - E/K) 2 ) p r 



(5.1) 



which is easily solved, putting [...] = const in eq. (5.1). Inserting the result in eq. ( |3 . 1 3^ we obtain the prefactors of 
the semiclassical wave function to lowest order. The exponential factors exp(±iS/h) in ( |3.l[ ) are combined to satisfy 
the usual matching conditions on the lower turning point. After some rearrangement in the prefactors we obtain 



u nr i m {r) 

Vn r lm{r) 



Cof VVE 2 + if 2 + K + y /yl 2 + K 2 - K 
2r ^ VV-E 2 +K 2 + K - VVE 2 + K 2 - K 



(5.2) 



where p r {r), K 2 (r), are defined in eqs. (4J), ( 4.4), an d r< = p<doi?o given by eq. (4^) for type B motion or ( 4.11 ) 
for type A motion. The normalization condition (2.13) serves to determine the constant Co by the integral 



1 = 



C 2 



dr 



E 



r 2 



VE 2 + K 2 \p r (r)\ 4m 



(5.3) 



where T r is the period of the radial motion, and r> = p>doi?o is defined by eq. ( 4.15| ). To obtain eq. (5.2) we 
have restricted the normalization integral to the classically allowed domain and r eplaced the rapidly oscillating 
sin 2 (?i~ 1 f r p r dr + n/4) by its average 1/2. Furthermore we used the identities ( |3. 18 ) for the radial velocity. 



A. General case 



The integral in eq. (5.3) can be carried out explicitely. The necessary steps are summarized in Appendix B. We 
obtain 



4mup 



c 2 



type A motion 
- i,. „,■_>. --.-i;*)-!/* type B motion 



4muo 



(5.4) 



with 
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a = arctan ■ 



(3 = arctan 



Ve - j 2 e + e 2 + 2j 2 + e-y/j + e 2 + 2fj 
v/2(e 2 + 2j 2 ) e 2 + j 2 + e^/l + e 2 + 2j 2 



(5.5) 



Using these results it can be shown that for type B-motion the wave-functions for large energy tend to those of the 
free harmonic trap. 



B. Large condensate 



The result (5.4) simplifies if we consider the limit Rq — > oo, keeping 



E= l -hiu Rle, £+1/2 = R 2 j 

fixed. Type A motion disappears in that limit and for type B motion we obtain 

2mo>o 



(5.6) 



Cl = f^J2 + h 2 cj 2 (£+l/2) 2 /E 2 

IT v 



(5.7) 

The correctly normalized wave functions ( |5.2| ) can now be used to evaluate matrix elements in the semiclassical limit. 
This will be done in the next section to calculate the cross section for inelastic light scattering by which a quasi-particle 
is created from the condensate. 



VI. SEMICLASSICAL STRUCTURE FUNCTION 



To introduce the dynamical structure function let us first consider the differential cross-section for the off-resonant 
scattering of light from the trapped gas. Light with incoming wave- vector and frequency q^, lj^ is scattered into an 
outgoing field with q' L = qi, — k,w' L — lul — oj with ui « lol- The scattering angle between q' L and q^ is related 
to k by the usual kinematical relation ck = 2iOi,sin(9 /2) which implies the allowed interval < k < 2ul/c. The 
permitted interval for the transferred frequency to can be estimated from the corresponding relation for free atoms 
< lo < (2TiLu\/mc 2 ). Let now the differential cross-section for the off- resonant light scattering from a single atom 
in its ground state be given by (da / d£L) At . Then the spectral distribution of the differential cross-section for light 
scattering off the trapped atomic gas is given by |l8||ll| 



d^ldOTiuj) \dfl 



(6.1) 



At 



The dynamical structure function S(k,co) is the Fourier transform of the density — density correlation function 
and defined as the thermal expectation value 



S(k, w) = ~X; e-^\(fj,'\p(k)\n) \ 2 S(hcu + E„ - E„) . 



Here \ and E^ are energy eigenstates and levels of the system, and 



p(k) = J d 3 rii + (r)e~ 



ik-r 



ip(r) 



is the operator of the density fluctuation with wave vector k. 
In the present case S(k, oj) assumes the form 



S(k, u) = S (k)S(fuj) + ^ ' S t {k) [S(hu + E t ) + e- fiE *S{Tvu) - E z )] 

i 

+53'4 1) ( fe ) [ S( - hiJ + E * + E i) + e-ft Ei+E >)5(hu -Ei - E 3 ) 



(6.2) 



(6.3) 



S§'(k)Q(Ei - Ej) 8{huo + Ei- Ej) + e-^-^Sfhu - Ei + Ej) 



(6.4) 
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where the sums J2' run over the quasi-particle states with Ei ^ 0. The first term in eq. (6.4) describes coherent 
clastic scattering, and So(k) is simply the square of the Fourier transform of the density 



d 3 re 



— ik- 



\vi(r)\ 2 +J2'(\* 



2 )C 



(6.5) 



Because tpo cx v^o this is by far the dominant contribution in S(k,ui). This contribution of the scattered light, by 
its interference with the non-scattered light, is used in the phase-contrast imaging of the condensates [^Tj. The part 
Yl'i \ v i( r )\ 2 gi yes the number dens ity of particles outside the condensate at temperature T = due to the interaction. 
It can be easily seen from eq. (5.2) that |ui(r)| 2 vanishes for r > tq. For weak interaction this number density is much 
smaller than Nq/tq. 



The second term in eq. (6.4) describes the creation of a quasi-particle from the condensate. The matrix element for 
this process is 



M*(fc) 



d i r e l "' r ip {r)(u*(r) - v*(r)) 



and 



(6.6) 



Si{k) = \Mi{k)\ 2 {ni + 1), n t 



The other quantities in (6.4) are: 



M#>*(fe) = - 



d 3 r ujvje 



ikr 



M l?*( k ) = I d A r{u* Ul +v*v 3 )e 



ikr 



Mg\k) y / (n l + l)(^ + l), sg\k) 



Dm 



Here is proportional to the cross-section for the absorption of an energy fuv by creating a pair of quasi-particles 
with energies Ei,Ej from the condensate. The cross-section for the time-reversed process is smaller by a factor 
e -/3(B i +_E j ) j ue to detailed balance. Similarly S^Q(Ei — Ej) describes the cross-section for absorbing an energy Tiuj 

by converting a quasi-particle of energy Ej into one with energy Ei. Both and give rise to a faint and broad 
spectral background in the scattering. In the following we shall therefore content ourselves with the evaluation of 
|M,(k)| 2 . 



(6.7) 



(6.8) 



A. General case 



Here we shall evaluate the matrix element (S.6) in the semiclassical limit. It is clear that only type B motio n can 
contribute, because otherwise ipo and (iii — Vj) don't overlap. The condensate wave function is given by eq. (2.5). The 
amplitude (it; — «,-) is obtained from eq. fl5.2| ), and e~ lk v , where k in z-direction in semiclassical approximation is 
expanded as 



lfe2 = £(- 



e 2 /7r(2£ + l) 



k 



k 2 




1/2 



y ro (e,$). 



(6.9) 



The int egra l ( |6.6[ ) is to be taken over the angles O, $ selecting the angular momentum quantum number from the 
sum in (3.5) and over r 2 dr from r = r < to r = tq. This radial integral is of the general form 



dr f(r) cos ip(r) 



2tt 



\<P"(*o)\ 



f(zo) cos (v(z ) + sgn (<p"{zo)) 



(6.10) 
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where the evaluation is done in the stationary phase approximation and zq is defined as solution of tp'(zo) = in the 
interval r< < zq < tq. Before writing the integral explicitely we introduce scaled parameters and variables as before 
and in addition 



k = Rqk/ da , Co = y mwgC 



(6.11) 



Then, neglecting a rapidly oscillating term without stationary phase in the interval r < < r < tq, the amplitude can 
be written as 



M nr i m (k) = S m0 (-iY\l ——2 



e JC 1_J15N (2£ + 1) f 1 



8 



max(j/K,p<) 



dpJl 



l-p 2 



Ve 2 + (1 - P 2 ) 2 



1 



: COS 



R? 



with 7T r (p) given by eq. (4.7) and 



dftir r {ff) ~ Rl / dp'wM 



The stationary point zq is determined by 



which is solved by 



k 4 + 2k 2 



2k 2 



Now we have to examine under which conditions zq lies in the required interval: 
i) The condition for ipo, u — v, and e~ lkz to overlap are 



Of 



J < Hi 



ii) zq satisfies < z 2 < 1 if 

iii) zq satisfies the lower bound j /k < zq if 



k 4 + 2k 2 > e 2 > k 4 



e 2 < k 4 + 2k 2 - 2]' 



(6.12) 

(6.13) 
(6.14) 

(6.15) 

(6.16) 
(6.17) 
(6.18) 



which is stronger than the first inequality (6.17) if the second inequality (3.17) holds, 
iv) zq satisfies the lower bound p < < zq if 

, ^ ^ e 2 (n 4 + 2Hi 2 -2j 2 -€ 2 ) + 2j 2 K, 4 , . 

y/l + e 2 + 2j 2 < 1 + -i —L- 1 J - . (6.19) 

ZJ K 

It follows from (6.18) that the right hand side of (3.19) is positive. Therefore both sides of ( |6.19 ) can be squared 
without changing the inequality, and we obtain after some rearrangements instead of ( 6. 19] ) 



f + 



2k 2 



2k 2 



< 



2k 2 



e 2 + K 4 
2k 2 



(6.20) 



The terms inside the brackets [. . .] on both sides of (6.20) are positive, by virtue of eq. (6.17). Therefore we ma y tak e 
the square root of both sides of eq. ( |6.20| ), to obtain after some more r earra ngements , again the inequality ( 3.18 ). 
Thus, eq. ( 3.19| ) leads to no additional condition if inequalities (3.16)— ( 6.18) ) are satisfied, which are therefore the 
semiclassical selection rules for the creation of quasi-particles from the condensate. For later reference it is usefull to 
rewrite the most restrictive of these inequalities in dimensional quantities as follows 
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2l2 



2m 



2m 



% 2 k 2 1 



m 



(6.21) 



A numerical example for this selsction rule is shown in Fig. |^. 
Next we need to determine tt w (zq) = ir r (z ) and obtain 



K 4 + 2K 2 -e 2 -2j 2 
k 4 + 2k 2 - e 2 



We also need the sign of 



(p (z ) = R (n' r (z ) - n w {z )) = R 



2k 2 



k 4 + 2k 2 - e 2 



e 2 + ^(k 4 + 2k 2 - e 2 - 2j 2 ) 



(6.22) 



(6.23) 



which is evaluated to be positive. We can now apply the formula ([11]) to ( |6 .11 ) and obtain 

M nrlm (k) = (~i) e 5 m0 



'C e /l57riV (2£+ 1)V2 



K K 3 R% 



\J k 4 + 2k 2 



2f 



: COS 



(^+4 



where the phase tp 

is given by the integrals 



(6.24) 



(6.25) 



</?i = Ro / dpir r (p) 
p-2 = -Rq / dpw w (p) . 



(6.26) 



Integrals of this kind had already to be evaluated when c alcul ating t he a ction I r , and the same substitutions as were 
made there are also of help here. In fact the integrals ( |6.26|) and ( 4.36 ) differ only by a pre facto r h/n and by the 
upper bou ndary . Using this observ ation we obtain by the series of steps summarized for eq. ( |4.36| ) in appendix A a 
result for ( 3.26 ) which differs from (A5) only by the prefactor and the boundary A\ — > A^ with 



A* = arcsin 



e - k 2 + e 2 + ek 2 + 2f 



(e + K 2 )^l + e 2 + 2y 



(6.27) 



Explicitely 



?! = Rli^2c i J > (A 3 ,a i ,bi) + C 4 J<(A 3 ,a4,&i) 



+C 5 K < (A 3 ,a i ,b 1 ) 



with the parameters given by eqs. ( |4.42| )-( |4.5l| ). The second integral ( p.26[ ) is simplified by the substitution 



KCOSU 



and gives, after evaluation 

^2 = -RqJ 



/k 4 + 2k 2 -£ 2 



2j 2 



1 — arccos 



V2j 



Vk 4 + 2k 2 - e 2 



(6.28) 



(6.29) 



(6.30) 



Now the matrix el eme nt is completely determined in the semiclassical approximation. The inelastic part of the 
structure function (3.4) at T — can be given by the sum 
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S (1 \k,u) 



irdnU) 2 



^/2a(kdo) 7 Ti(jj, 



C nr£ (2l+l)[(2Lj/u ) 2 -(kd ) 4 



(2uj/u ) 2 -2(£- 



(kd ) 4 
1 — sin 2(p 



2{kd Q fRl 



5 



W + U Hr i 
W 



(6.31) 



where the two sums are restricted to quantum numbers n r , £ satisfying the semiclassical selection rules ( 6.16|)-(|6.18] ) 
with e = huj nr g/[j, and j = (£ + ^)/Rq. Here we have put indices n r , £ on C and <p defined by eqs. ( 6. 11 ), ( [5. 4| ) and 
(6.25), respectively, to indicate the dependence of these quantities on the quantum numbers n r ,£ via e = hu n g/fj, and 
j = (£+±)/R 2 . 



B. Limit of large condensate 



In the limit of a large condensate R Q » 1, E/huj a ~ 0(1), (£+1/2) ~ 0(1), kR d ~ 0(1) the result ( |6.31| ) f or th e 
inelastic structure function simplifies. For the prefactor C Ur i we can use the asymptotic form (5.7) with eq. (3.11). 
The semiclassical selection rules (6.16)- fl6.18| ) take the form 



1 (^+l/2) 2 

E nr e > -jTlUJQ ^ 



1/2 < R a kd 



(6.32) 



1 



1 



^o(fcdo) 2 < E nrt < ^=hio QX / Rl(kdo) 2 -(£ + 1/2) 2 



(6.33) 



where the right-hand sid e of (16.33 ) has been simplified using (kdo/Ro) 2 << 1. 

The first condition of ( |6.32| ) is always satisfied for sufficiently large Rq, the second represents a restriction on the 
size of £ since we want to assume Rokdo ~ 0(1). In (6.33) we can use the result (4.25) for the quasi-particle energies 
in large condensates and obtain 



2n r +£ + 3/2 < R kd 



(6.34) 



for the second inequality in (3.33), while the first is al ways s atisfied in the limit we consider. It is somewhat surp rising 
to see the energy levels of the free trap appearing in ( |6.34 ). Under the square-root in the denominator of eq. (6.31) 
and als o in t he nom inator we may ag ain neglect (fcc?o) 4 compared to the other terms. Finally, the tedious general 



results 



(4.18), (4.19) for p< and 7r r , and 



.28), ( 5.30 ) with ( 4.40 )-( 4.51 ) for tpi and ip% may be simplified considerably by using the asymptotic results 



z -> yl - 2(uj n jJuoRokdoY 



(6.35) 



MP) = ^L^(R a kdo) 2 -{£+l/2f/p 2 



R l 



(6.36) 



following from eqs. ( |6.15 ), (|6.13 ) in eqs. ( |6.26| ). The integrals are now much easier to perform and we obtain after 
some further rearrangements 



<Pn r t(k) = -Rokdo (sinxn r e(k) - Xn r e{k) cosxn r e(k)) 
with the angle Xn r i(k) defined by 

2n r +£ + 3/2 



cos Xn r e{k) 



Rokdo 



(6.37) 



(6.38) 



The inequality (6.34) is incorporated in this definition. After all this the structure function QS.31 ), in the limit of 
large condensates, takes the explicit form 



5 (1) (fc,w) 



d 



afioJo (kdo) 7 u)q 



(2n r +£ + 3/2)(2£+l) 
„ , ^ (Rokdo) 2 - (2n r + £ + 3/2)= 

(2n r +£+3/2<fi fcdo) 



EE 



(6.39) 



cos" 



+ Vn r i(kfj 5 ((uj + Lo nr£ )/ujo) 
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with ip nr i and ui nT i given by eqs. ( |6.37 ) and ( 4.25| ) respectively. It is interesting to note that SW(fc,w) in the present 



limit has resonances for the maximum momentum transfer Rokdo = 2n r + £ + 3/2 by which states with the principal 
quantum number n = 2n r + £ ca n be excited. The angles Xn r e(k) and ip nr e(k) vanish there. Similar resonance s app ear 



also in the more general result ( 6.31 ) which diverges whenever the border of the semiclassical selection rule ( 6.1S ) is 
approached. 

VII. SOME NUMERICAL EXAMPLES 

All the calculations, so far, have been analytical, but the final results for the energy levels, wave functions and 
inelastic structure function are not obtained in explicit form and must be evaluated numerically from the implicit 
results we have derived. This can be done without difficulty. We shall here consider some examples of such results: 

A. Energy levels: 



The energy levels are obtained by solving eq. ) for E. This can be done numerically without difficulty and 
ambiguity, since 



dl r _ T r 
~dE ~ 2^ 



>0 (7.1) 



where T r is the period of the radial motion. In Fig. || we take (/i/Hujo) = ^R 2 = 4.27 and plot the level shifts 8 n (£) 
defined as a function of £ and the principal quantum number n = £ + 2n r by 

huo5 n (£) = E nrt + (i - hcj (n + 3/2) . (7.2) 

As for type A motion the levels are the sam e as for the free trap, apart from a shift in the zero-point energy from 



3hujo/2 to 3hujo/2 — /x, it is clear that (7.2) vanishes for this type of motion, characterized by the simultaneous 
inequalities 1 < 2j — 1 < e < j , which imply 

{E nri + n)/hu >(£+ 1/2) > {2E nrl /hoj ) 1 ' 2 {2^/hu Jo ) 1 ' 2 . (7.3) 

The first inequality must al ways be satisfied for energy levels to exist, at least semiclassically, i.e. there are no levels 
where it is violated (cf. eq. ( 4 . 1 2| ) ) . The second inequality marks the border in the (n, ^)-plane in Fig. || where 8 n {£) 



drops to 0. For type B motion £ r + 1/2 < (2E rir e / Kuio) 1 I 2 (2/j,/ Hujo) 1 ^ 2 is satisfied. The level shifts are positive due to 
the repulsion provided by the condensate. This repulsion is stronger for the lower lying levels, as one could expect. 

B. Structure function: 

The dynamical structure function in the approximation we have considered is determined, up to prefactors, by 
just two parameters, Rq which fixes the number of particles of a given atomic species in a given trap via Nq — 
\J % j 'mtuoRo I '15a, and kdo which fixes the scattering angle via 9 = 2Arcsin [fcc?o(A/47r)(mwo/?i) 1 / 2 ] , where A is the 
wavelength of the scattered light. We shall assume that the wavelength of the light used in the scattering experiment 
is off-resonance but roughly given by the D-line of the alkali-metal, i.e. A ~ 590 nm for sodium, and A ~ 800 nm for 
rubidium. 

The result ( |6.3l| ) for the inelastic part of the structure function for R$/2 = fi/huiQ — 4.27 and kdo — 15-1 is plotted 
in Fig. ||. For a trap with frequency vq = luq/2tt = 100 Hz these numbers correspond to light scattering from a 23 Na 
condensate of N Q w 8000 atoms with 9 » 40° or from a 87 Rb condensate of 7V « 4000 atoms with 9 w 126°. This 
is a typical example for light scattering with a large momentum transfer in which high -lying quasi-particle states are 



excited. The plot gives as discrete points the strength of the 5-functions in eq. ( 6.3l| ), which according to our sign 
convention for u> appear at negative values to = —E rir e/h for absorption of energy. Since we consider temperature 
T = in the plots of Figs. [| the emission lines at positive u = E nr e/fi are frozen out. The infinitely sharp lines 
described by the 5-functions are, of course, not directly observable. Apart from the intrinsic line-widths which we 
did not consider in this work, they must also be convoluted with the experimental resolution. However, in order to 



retain the full information present in eq. ( 5.31 ), we prefer not to carry out such an arbitrary smoothing. The inelastic 



part of the structure function is seen to consist of several discrete quasi-continuous rotational bands, roughly spaced 



17 



by ?uvo. In the region of parameter space chosen here, these bands can be usefully labelled by the principal quantum 
number n = 2n r + £ of the free trap, in terms of which the energy levels are expressed as 



E nrt = huj {n + 3/2) -fi+ huj S n {£) . (7.4) 



If 5 n {£) = were satisfied each band would collapse to a single sharp line. The semiclassical selection rule (6.17) 
limits the excited bands to the domain 

J/W^Hf >f^>^. (7.5) 

In the present case (kdo) 2 >> 2Rq we may simplify this relation somewhat by expanding the square root and using 
Rq = 2li/TujJq to obtain an energy interval of size li for the excited quasi-particle states which extends above the recoil 
energy h 2 k 2 /2m 

(kd ) 2 + _V_ > ^nA > (kd ) 2 
2 huiQ ~ hu! ~ 2 

In the local density approximation (22) the structure function is just a single continuous band in the same energy 
interval (and a corresponding band at negative energies related by detailed balance). It is not surprising that our 
treatment, being based on the discrete quasi-particle sta tes, reveals considerably more structure, namely the bands 
and their discrete substructure. As can be seen from (7J3) the number of bands in Fig. § just gives the integer part of 



ll/Tiluq. This rather direct physical measure of li in units of hu>o could be of practical ex perim e ntal i nterest. The fine 
structure of the rotational bands is also determined by the semiclassical selection rules ( |6 . 1 6| ) ( 6.18 ). The rotational 



quantum number £ increases within each band of Fig. || from le ft to right. The strongest line within each band 
corresponds to the maximum value £ m ax(n) of I permitted by eq. ( 6.1 8| ) for type B motion. 



We can compare these results with similar ones one may infer from Fig. 2 of ref. |14| | , w here the quasi-continuous 
rotational bands due to high-lying quasi-particle states were first found. The results in |H] were obtained by a much 
simpler approximation, not aiming at achieving full self-consistency, using the excited state wave functions of the 
empty trap, rather than the WKB wave functions whose use underlies the present Fig. |[ This implies that ir r (p) in 



eq. (3.12) is replaced by the free oscillator expression 

4 = Vt-.f/P 2 -P 2 - (7-7) 

Of course this leads to much simpler closed formulas for |M n ^(fc)| 2 (eqs. (5), (6) of ref. [Q). However, with our 
present self-consistent treatment we are now able to assess the limitations of the simpler approximation. In fact, the 
result presented in Fig. || differs in two important details from that of jTJj : 

(i) There is an overall shift of size /i of the energy window excited by the transferred wave number k, and (ii) there is 
a qualitative change in the dependence of the width of the rotational bands on the principal quantum number which, 
according to Fig. |3J, become more narrow as the principal quantum number increases, while the oppo site b e havio r 
was found in JlJ]. Both differences have their origin in a change in the semiclassical selection rules ( j6.16| )— ( |6.18 ) 



which depend rather sensitively on approximations in the phases of the wave functions inside the condensate. The 
first shortcoming of the approximation in [|l4j can be easily eliminated by using simply E° sc — fiu Q {n + 3/2) — fx 
rather than the free trap energies in the expressions for the free trap eigenfunctions. The second is not as easily 
avoided, however. It simply shows that the influence of the interaction with the condensate on the spatial dependence 
of the phase of the high lying states cannot be neglected as long as the widths A„ of the rotational bands given by 
A„ = <5„(0) - S„(£ max (n)) with 

£ m M + 1/2 = V2 ((fcd ) 4 /4 + (k<k) 2 Li/huj - (K s 7^o) 2 ) 1/2 (7.8) 

are not themselves negligible. We conclude that for a proper treatment of the rotational bands the present treatment 
is indispensable. 

The spectroscopic resolution of the quasi-continuous bands requires a resolution better than the widths A„ , which 
may be difficult to attain in the frequency domain, but could be feasible in real-time experiments extended over a 
time-scale of w 10 periods 2it/ljq. In real time the width A„ would turn up as collapse rate 7„ = A„ determining the 
dephasing of the band with principal quantum number n. Due to the discrete substructure of the bands this dephasing 
should in principle be reversible. However it appears unlikely that the strict phase coherence over a time-interval 
of the order of £ max (n) / \ n , required to see such a rephasing, could be achieved experimentally. For a discussion of 
collapses and revivals of low-lying collective modes see E§] |29] . 
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Let us now look at the dynamical structure function also for smaller momentum transfer and for somewhat larger 
condensates. In fig. || we have chosen the parameters i? 2 , = 11.24, kdo — 3 corresponding to, No — 6000, 9 = 21° 
for &7 Rb and to No = 12000, 9 = 7.7° for 23 Na, if the same choices are made for vq and A as before. The changes 
in the scale in fig. ^ in comparison with fig. [| should be noted. In fig. ^ we have joined spectral lines with equal 
principal quantum number n — 2n r + £ by straight lines. This is uscfull because the quantum number n organizes 
the spectral lines according to multiplets of different strengths. However, differently to the case of high lying states, 
where n labelled well-separated quasi-continuous bands, the different multiplets now overlap at least partially in their 
energy range. The strongest multiplet in fig. ^ corresponds to n = 9, the next to n = 8, and the third and weakest 
one shown to n — 7. There are actua lly still weaker lines with n < 6 which we have suppressed. All the lines satisfy 
the semiclassical selection rule ( |3.21 ). It is plotted for this case in fig. ||| in the (angular momentum, energy)-plane, 
together with the semiclassical energy levels. All levels within the curved triangular section of the plane satisfy the 
selection rule. In this plot multiplets of fixed n = 2n r + £ lie on lines roughly parallel to the upper border of the 
triangular shape, the higher lying l ines corresponding to higher values of n. On this upper border the square-root in 
the denominator of the expression (6.31) for S^(k,u>) vanishes, leading to a scattering singularity, while on the lower 
border the nominator of this expression vanishes, leading to a suppression of scattering. This observation explains 
qualitatively the dependence of the scattering cross-section on n. The multiplet with n = 9 consists only of the 
two lines for £ = 1, £ = 3. The remaining £ = 5, £ = 7, £ — 9 members of this multiplet have smaller energy but 
lie above the triangular region in fig. 0. This highest lying multiplet of fig. H is therefore rather narrow and well 
separated in energy from the other multiplets. It lies entirely above the chemical potential /i, i.e. its lines in fig. || lie 
at oj/ujo < fi/Jiuio — —5.62. With these features this multiplet bears some resemblance to the well-separated bands 
of fig. |[ In contrast to this the two lower lying and weaker multiplets in fig. || overlap completely with each other 
and extend over energy levels ranging from values larger to values smaller than /i and already resemble the low lying 
multiplets extending down to the hydrodynamic regime. This can be seen from fig. [}] where we have chosen Rq = 26.2, 
kdo = 2 corresponding to N = 50000, 9 = 14° for 87 Rb and N = 10 5 , 9 = 5.2° for 23 iV a, again assuming a trap of 
vo = 100 Hz and off-resonance scattering at wavelengths roughly equal to the respective Z?-lines. The multiplets are 
still organized according to their strengths by the principal quantum number n. However, now these multiplets are 
completely overlapping, much as the two lower lying ones in fig. ^. For the multiplets shown in fig. ^| n ranges from 
n = 8 for the strongest multiplet to n = 6 for the weakest one shown in the figure, but again we have suppressed 
in the plot still weaker lines corresponding to multiplets with smaller values of n. The dependence of the scattering 
cross-section on n, and within a given multiplet on £, can be understood similarly to our discussion of fig. |J by using 
the asymptotic expression (6.38) for S^ 1 ' (fe, u>). Taken together the plots fig. ||-fig. ^ show that a measurement of 
S(k, u>) over a range of scattering angles 9 in principle can give a very detailed information on the quasiparticle levels 
in any energy interval selected by the choice of k according to eq. (6.21). 



VIII. CONCLUSIONS 



In the present paper we have solved the Bogoliubov equations for the elementary excitations of a Bose-condensed 
gas of atoms trapped in an isotropic harmonic potential in WKB-approximation. The underlying classical dynamics is 
integrable in this case due to the rotational symmetry. The quasi-particle energies in the trap have been obtained from 
a Bohr-Sommerfeld quantization rule. Corrections to that rule due to the non-differentiabilty of the condensate wave 
function in Thomas- Fermi approximation have been discussed in [^5|. The mode-functions, including the Bogoliubov 
amplitudes Vj , have been obtained in the classically allowed regime in the WKB-form. Our results interpolate between 
the low-energy regime E / fi — > 0, where asymptotically exact results within the Bogoliubov approximation are available 
H and the high-energy regime \ijE — > 0, where perturbation theory can be used. As an application we have given 
a self-consistent calculation of the fully resolved dynamical structure function S(k,ui), which could be measured in 
off-resonant light scattering, in a large energy domain extending to all but the lowest lying levels. 
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APPENDIX A: EVALUATION OF THE ACTION INTEGRAL 



Here we summarize the essential steps for the evaluation of the integral in eq. ( 4.36 ). The integral is splitted into 
two parts, the first part /< taken over p from p < to 1, the second part 7> taken over p from 1 to p > . Let us consider 
the evaluation of /< first. The substitution 



2et 



dp = —e 



(1+i 2 ) 



(l-t 2 ) 3 /2(i_ i 2_2ei) 1 /2 



dt 



transforms /< into an integral over t from 

t = to t = t raax 



-(e 2 +j 2 ) + e^l + ^ + 2f 



e + f 

A partial fraction expansion of the integrand results in the explicit form 

TiR 2 



I< = 



e 2 + p 



e + j" 



V^TT + e- 1 



(e + f) 2 
1 y/^TT-e+l 



1 



4e 2 
1 1 



V / e 2 Tl-e-i 
e- 2 1 1 



4e 2 



Ve 2 + 1 + e + t 



Ael-t 4e 2 1 + t 2e(l + tf 



A further substitution 



t = 



2f 



sm 

and treating the curly bracket in eq. (|A3|) term by term reduces the integral to the form 



e 2 +J 2 



I< = ^<j^C l J>(A 1 ,a l ,6 1 ) + C 4 J<(^i,a4,6i) 
,i=i 

+C 5 K < (A ll a i ,b 1 ) \ 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



with the functions J<, J>, K < defined in eq. ( 4. 37 )— ( 4. 38 ) and the parameters and coefficients given in eqs. ( 4.4C ) 

(EH). 

Turning to the second part i> of I r we may apply a similar sequence of substitutions resulting in 



P = 



e+1 



e + 1 



j z sin 99 



(A6) 



which brings I> to the form 



7> = ^lc 6 J > (A 2 ,a 5 ,b 2 ) 

7T 



with the parameters and coefficients given in eqs. (4.40)^(4.52). From I r =/<+/> we obtain the result given in 
eq. (|0|). 
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APPENDIX B: NORMALIZATION OF THE WKB WAVE FUNCTION 



Here we summarize the necessary ste ps i n the evaluation of the integral in eq. (5.3). For type A motion t he radial 
period is simply tt/ujq. To evaluate eq. (5.3) also for type B motion we introduce scaled variables like in eq. (4-6), and 
the integral is splitted in a part I± from p < to 1 and a part I2 from 1 to p > . 



h = 



h = 



1 



dp 



1 



1 



dp 



mojQ J 1 

A series of substitutions amounting to 



e + 1 - £ - p 2 



P = 



1 - V - 2et 
1-t 2 

e\J\ + e 2 + 2j 2 sin^ - e 2 - j 2 



puts I\ into the form 



h 



I 1 ^ /•t/2 



dip 



mujQ J Al e + e 2 + 2j 2 - e^/l + e 2 + 2j 2 sin ip 



where yli is defined in (4.40). The integral is now elementary and given by 

2e /3 



h = 



v / 2(e 2 + 2f) muj Q 



with j3 given by eq. ( p.5[ ). Here we used the addition theorem for arctan to simplify the final expression. 
I2 is much simpler to evaluate. The substitution 



reduces it to the form 



2 e+1 

P = ~^r- 



h = 



e+1 



j 2 sinyi 



2mujQ 



it/2 



dep 



A-2 



where Ai is defined in eq. (4.41). After some simplification the result takes the form 



h = 



1 /7T 



2mujQ V 2 



(Bl) 
(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



with a given by eq. (5.5). Combining the results flB5|), (B8) with eq. (5.2) we obtain eq. (p.4|). 
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FIG. 1. Regions of type A motion and type B motion in the scaled angular momentum-energy (j — e)-plane. In region D 
there is no lower turning point of the radial motion, i.e., in the semiclassical quantization region D is forbidden. 

FIG. 2. Difference S — S n (£) of quasi-particle energy levels, in units of htuo, to the levels of the free trap. The chemical 
potential was chosen to be fi = 4.27fttJo. 

FIG. 3. Dynamic structure function S^(fc,aj) in units of do/ahuio for fi/hujo = 4.27, kdo = 15.1 as a function of to in units 
of u>o- Plotted are the strengths of the Dirac-delta peaks. 

FIG. 4. Dynamic structure function S^ (k, lu) in units of do/ahuio for fi/hojQ = 5.62, kdo = 3 as a function of ui in units of 
ujo- Plotted are the strengths of the Dirac-delta peaks. Peaks with equal principal quantum number n — 2n r + £ are connected 
by straight lines, with n = 9, 8, 7 from top to bottom. 



FIG. 5. Semiclassical selection rule 6.21 and semiclassical energy levels E — E n ,,i in the (£, L?)-plane for fi/huo = 5.62, 
kdo = 3. 
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FIG. 6. Dynamic structure function S 1 ' 1 ' (k, oj) in units of do/ahuio for h/Uojo = 13.1, kdo = 2 as a function of uj in units of 
uio. Plotted are the strengths of the Dirac-delta peaks. Peaks with equal principal quantum number n = 2n r +£ are connected 
by straight lines, with n — 8, 7, 6 from top to bottom. 
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